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ABSTRACT 

The  log'X^  distribution  of  correlation  times  and  its  application  to  polymer  systems  are 
examined.  It  is  found  that  the  original  mean  correlation-time  parameter,  xq,  (Macromolecules, 
1973,  6,  882)  does  not  represent  the  true  mean  of  the  distribution.  A  properly  parameterized  log- 
distribution  and  function  relating  the  original  Xq  and  width  parameters,  p  and  b,  to  the  true  mean 
are  derived.  Fitting  of  both  the  original  log-x^,  corrected  log-x^,  and  log-normal  distribution  to 
experimental  data  is  discussed  and  the  log-normal  distribution  being  proposed  as  an  alternative. 
Interpretations  of  NMR  data  from  polymer/solvent  systems  can  be  different  for  the  log-X^  and  log¬ 
normal  distributions.  Due  to  its  better  defined  moments,  the  log-normal's  results  are  believed  to  be 
more  physically  realistic  when  large  widths  are  required. 
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Often,  definitive  models  for  molecular  motion  do  not  exist  or  fail  to  give  reasonable  results. 
And,  as  an  approximation,  an  exponential  decay  of  correlation  with  a  distribution  of  coirelation 
times  is  used  as  the  correlation  function.  Such  distributions  prove  very  useful,  especially  for 
polymers  and  polymer  solutions,  at  describing  and  quantifying  the  dynamics. 

Heatley^  has  summarized  the  successes  and  failures  of  several  of  the  distributions  used  to 
describe  polymer  dynamics.  Of  these,  the  asymmetric  log-x^  distribution  has  had  good  success  in 
fitting  data  from  a  bnoad  range  of  techniques  including  NMR,  ESR,  and  dielectric  studies.  ^ -2.5.8 
We  have  found  that  the  use  of  the  log-x^  distribution  may  occassionally  yield  parameters  which 
should  not  be  compared  from  system  to  system.  In  some  cases,  a  proper  comparison  can  be  made 
with  the  application  of  a  correction  factor.  For  example,  the  application  of  the  log-x^  distribution, 
as  proposed  in  the  literature,^  to  NMR  relaxation  data  of  toluene  in  polystyrene  generated  a  slower 
mean  rotational  correlation  time  for  pure  toluene  than  it  did  for  motionally  hindered  toluene  near 
polystyrene  at  the  same  temperature.  It  also  yielded  unphysical  results  in  a  dynamics  study  of 
polyO'jo-propyl  acrylate)  (PIPA).2  In  that  study,  the  mean  correlation  time  predicted  by  the  log- 
X^  distribution  had  little  correlation  with  polymer  concentration,  often  being  faster  for  more 
concentrated  solutions  than  for  less  concentrated  ones  at  the  same  temp^ature.  Other  models 
employed  to  quantify  the  PIPA  system  suffered  similar  problems.  After  examining  the  log-x^ 
distribution,  we  found  that  the  parameter  Tq,  previously  inteipreted  as  the  mean  correlation  time,  is 
not  the  true  mean,  but  rather  a  scaling  factor  which  can  be  related  to  the  true  mean. 

A  normalized,  scaled,  base  b  log-x^  distribution,  FX(Xc;  b,p,0),  has  the  form  of: 


FX(x,;b,p,e)  =  l!>M 
e  r(p) 


for  O^Xc^oo 


(1) 


with  Ip,  b,  p,  0,  and  F  the  correlation  time,  numerical  base,  width  parameter,  scaling  parameter, 
and  gamma  function,  respectively.  The  superscript  x  denotes  the  log-x^  distribution.  The 
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panuneter  b  is  often  fixed  at  1000,  but  it  could  be  replaced  with  e,  the  natural  base.  The  mean, 
'^(b,p,e),  of  FX(Xc;  b,p,0)  is; 


too 

).p.8)=j  xcF«(tc:b,p,e)dii;  =  e[jp^]f’-i] 


Experiments  often  probe  the  dynamics  through  the  Fourier  transform  of  the  correlation 
function  or  spectral  density,  J(o)).  For  a  distribution,  F(X{.),  of  exponentially  decaying  correlation 
times,  J((D)  is: 


Jfoo 

I 

0  1  + 


aP-xi 


where  <0  and  Xg  are  the  frequency  and  correlation  time,  respectively.  FX(Xc;  b,p,6)  could  be  used  in 
equation  (3)  and  the*  mean  correlation  time  calculated  from  equation  (2)  after  determining  the 
parameters  b,  p,  and  6.  However,  transforming  the  distribution  and  spectral  density  to  a  reduced 
coneladon  time  will  allow  a  more  direct  comparison  with  the  original  work^  and  eliminate  the  need 
for  equation  (2)  as  the  scaling  parameter,  6,  will  be  replaced  by  a  mean  correlation-time  parameter. 


With  the  transformation  x^  =  Xc/^(b,p,6)  applied  to  FX(Xc;  b,p,0),  the  log-x^  distribution  of 
reduced  correlation  times,  G^^fXr;  b,p),  is: 


_(ln(b)r 


GZ(Xr;  b,p)  =y2^tj(b,p,l)  (ln(XrtJ(b,p,l)  +  1)^* U^fb.p,!)  +  (4) 

r(p) 

For  simplicity,  0  is  set  to  one  in  the  arguement  as  the  transformation  eliminates  any  0  dependence. 
Similarly,  J(co)  becomes: 
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«.  f 

Jo 


2GX(Tf;b.p)TrTo 
1  +£l)^X?to 


(5) 


where  Xq  has  rqplaced  T^(b,p,6)  and,  consequently,  represents  the  true  mean  of  the  distribution. 

If  FX(xc;  b,p,6)  had  been  scaled  by  the  parameter  Tq,  as  in  reference  (1),  instead  of 
x§(b,p,0).  To  of  equation  (5)  would  not  be  the  true  mean,  but  a  scaling  parameter.  The  mean 
correlation  time  could  then  be  calculated  from: 


Xc  FX(Xc;  b,p,0)  dXc  =  Xq  j  x,  G'X(Xr;  Xo,b,p,0)  dXr  (6) 

Jo 

where  G'Z(Xr;  Xo,b,p,0)  is  FZ(Xc;  b,p,0)  scaled  about  the  parameter  xq.  Being  scaled  about 
x§(b,p,0)  simplifies  matters  by  cancelling  any  previous  scaling  and  making  the  mean  of 
G5C(Xr;  b,p)  =  1  and,  so,  <Xc>  =  Xq. 

For  comparison,  the  spectral  density  (equation  (5))  can  be  but  into  a  form  similar  to  equation 
(25)  of  reference  (1): 


J(o);  x^b.p)  = 


2x^(s,P)[eS-l]xo 


ds 


(7) 


where  p,  by  and  09  are  as  defined  previously,  s  is  an  integration  variable, 

X^(s,p)  =  -2L  sP-1  e-s  P  (8) 

r(p) 

and  Xq  is  the  mean  correlation  time  of  the  distribution. 

The  relationship  between  the  true  mean,  xq,  of  equations  (5)  and  (7)  and  the  originally 
reported  mean,^  can  be  calculated  from  equation  (6)  and  is: 
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pt^g  j 

^  =  -  (9) 

where  the  superscript  org  indicates  parameters  that,  were  fitted  from  the  original  distribution. 
Depending  on  b  and  p,  this  could  be  a  large  correction  and,  consequently,  responsible  for  the 
log-x^'s  previously  reported  unphysical  behavior.  Also,  this  equation  and  the  distribution  itself 
are  only  valid  for  p  >  ln(b)  where  the  distribution  is  normalizable.  The  possible  limitations  of  the 
distribution  and  its  dependencies  on  p  and  b  are  clarified  by  examining  the  moments  of 
GX(Tr;  b,p)  which  are: 


(10) 


where  ( ?)  is  the  binomial  coefficient  The  original  distribution  has  a  similar  form  for  its  moments. 
For  large  values  of  p,  the  moments  are  essentially  independent  of  b.  For  the  intermediate  values, 
ln(b)  scales  p  to  smaller  values,  making  the  distribution  wider  for  a  given  p/ln(b)  ratio.  Finally, 
for  small  values  of  p,  the  dependence  is  unclear  as  the  moments  are  infinite.  This  behavior 
suggests  that  the  natural  base  e  should  work  as  well  as  any  other,  i.e.  1000,  as  long  as  p  is  not 
restricted  to  integer  values.  Indeed,  with  non-integer  p  values,  all  data  tested  fit  equally  well  and 
with  the  same  predicted  Xq  using  either  e  or  10(X)  as  the  base. 

To  test  the  corrected  log-x^  distribution’s  ability  to  fit  polymer  data  and  yield  physically 
meaningful  results,  we  used  it  to  fit  four  data  sets:  polystyrene,  polyisoprene,  and  polybutadiene 
T1-T2  (20.6  MHz)  data  from  the  original  paper^;  polyisoprene  Ti  (20.1  and  100.6  MHz) 
data^;  poly(«o-propyl  acrylate)-d  in  chloroform  ^H  T1-T2  (13.8  MHz)7  data;  and,  toluene-dj  in 
polystyrene  ^H  TpT2  (30.7  MHz)  data.^  These  systems  represent  a  wide  variety  of  polymer 
dynamics,  namely  bulk  polymer  dynamics,  solution  state  polymer  dynamics,  and  solvent  dynamics 
in  polymer  solutions. 
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For  the  polymer  data  in  the  original  paper,  the  mean  correlation  times,  Xq's,  from  both  the 
original  and  modified  log-x^  seem  physically  reasonable  for  polystyrene  in  o-dichlorobenzene  and 
bulk  polyisobutene.  However,  the  original  distribution  yielded  an  unusually  fast  Xq  of  0.1  ns  for 
bulk  polybutadiene.  Using  equation  (9),  this  ‘q)  became  5  ns,  the  same  value  Ot  from  the  modified 
distribution.  This  seems  more  reasonable.  In  general,  the  modified  log-x^  distribution  generates 
larger  Xq's  than  the  original.  But,  with  the  use  of  equation  (9)  both  distributions  yield  the  same 
Xq,  with  smaller  values  of  p  requiring  larger  corrections.  These  results  are  typical  of  all  the  data 
fit  in  this  paper. 

Unlike  T1-T2  data  investigated,  both  the  modified  and  original  (in  conjuction  with  equation 
(9))  log-x^  distributions  generate  unrealistically  long  Xq's  from  the  polyisoprene  Tj  data.^ 
Table  I  lists  the  results  from  the  data.  The  reason  the  log-x^  distribution  fails  for  the  data 
is  not  clear.  It  may  linked  to  its  undefined  monwnts  and  that  the  data  only  samples  the 
distribution  in  the  MHz  region.  As  can  be  seen  from  Table  I,  the  modified  log-x^  compares  well 
until  p  becomes  less  than  2*ln(6)  below  which  only  the  log-x^'s  first  moment  exists.  For  these 
values  of  p,  the  distribution's  tail  is  poorly  parameterized  and  large  correlation  times  may  be 
unduely  weighted.  To  successfully  fit  data  to  such  a  distribution,  some  constraints  on  the  tail,  as 
the  zero  frequency  T2  data  imposes,  would  probably  help  give  more  reasonable  results. 

As  an  alternative  to  low  frequency  data,  the  log-normal  distribution  might  be  used.  Its 
moments  exist;  so,  its  tail  and  skew  are  well  defined.^  In  general,  it  has  simple  dependencies  on 
its  parameters  (easy  to  fit);  is  well  behaved  with  well  defined  moments;  and,  is  more  skewed  than 
the  log-x^.  It  also  affords  a  simple  physical  interpretation  of  the  underlying  processes. lO’U, 12 
The  distribution  of  reduced  correlation  times,  G(Xr;  a),  for  the  log-normal  distribution  is: 

G(v,  o)  =  1 —  e-(ln(Xr)/o-Ky/4)2  (1 1) 

Vic  axr 

which  would  replace  GX(Xr;  b,p)  in  equation  (5).  Its  second  moment  is: 


(12) 
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and  its  mean  is  one,  as  expected  firom  the  scaling. 

The  log-normal  distribution  fit  both  the  T1-T2  and  Tj  data  discussed  ealier  with 
vanishing  error  and  few  iterations  of  the  fitting  algorithm.  For  these  data,  except  for  where  log- 
failed,  the  Xqs  were  essentially  identical  between  the  modified  log-x^  and  the  log-normal, 
differing  from  the  original  log-x^  according  to  equation  (9).  As  mentioned  previously,  the 
predicted  magnitudes  and  trends  of  the  Xo's  ^cre  more  realistic  than  those  of  the  original  log-x^ 
and,  for  the  PIPA  data,  also  more  realistic  than  other  models  employed^  (compare  Figure  1  to 
similar  plots  in  reference  (7)).  That  is,  the  Xq's  increased  with  increasing  polymer  concentration  or 
decreasing  temperature,  realistic  in  an  intuitive  sense. 

For  reference  (5),  the  main  difference  between  the  log-normal  and  the  original  log-x^  is  that 
the  log-normal  predicts  a  higher  apparent  activation  energy  for  the  high  temperature  data  and  a  10 
fold  increase  in  Xq  at  low  temperatures.  As  in  the  reference  (5),  the  NOEF's  were  not  used  in  the 
fitting  procedure.  They  offer  an  independent  test  of  each  distribution’s  predictive  ablilty. 
Unfortunately,  the  calculated  NOEFs  follow  the  exnerimcntal  values  for  all  the  distributions  and 
no  distinction  can  be  made.  For  extremely  bro^  distributions,  however,  the  log-normal  has  a 
limiting,  Xq  independent,  NOEF  value  of  1.33  compared  to  the  log-x^’s  value  of  0.9  (at  20  Mhz). 
This  increase  is  presumably  due  to  the  larger  skew  of  the  log-normal  distribution.  None  of  the 
experimental  data  required  broad  enough  distributions  to  test  this  difference. 

In  conclusion,  Xq  original  log-X^  distribution  is  a  scaling  parameter  which  is  often,  but 
not  always,  representative  of  the  mean.  Comparisons  and  predicted  trends  may  be  on  different 
scales  unless  a  correction  to  the  true  mean  (equations  (6)  and  (9))  is  made.  Modifying  the 
distribution  yields  a  new  distribution  in  which  Xq  is  the  true  mean.  However,  both  the  original  and 
modified  log-x^  distribution  have  limited  use  as  their  moments  can  be  ill-defined  and  careful 
sampling  of  the  frequency  space  is  required.  We  prefer  the  log-normal  distribution  as  it  is  better 
skewed,  and  has  well  defined  moments  and  simple  derivatives.  Its  ability  to  fit  and  give  physically 
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reasonable  results  from  low  and  high  frequency  T1-T2  data  as  well  as  soley  high  frequency  T  j  data 
indicates  that  it  may  useful  where  the  log-x^  appears  unphysical.^«^*®**^*^5 
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Tabie  1 

Fits  to  the  Ci  methylene  Tj  data  of  reference  (5)  (Mn  =  7  kg/mol)  {b  =  1000) 

log-Xorg^  tog-Xwg*^  tog-Xmod^^  log-normal^ 

T  Ti(lOOMHz)  Ti(20MHz)  ToP  “lOP  “toP  XqO 


tsl 

fsl 

tns) 

mmu 

293 

0.152 

0.035 

2.00 

7 

1.90  7.3 

340*103 

7.5 

20.0 

3.45 

303 

0.156 

0.043 

1.00 

7 

0.67  8.0 

5500 

8.0 

7.10 

3.36 

313 

0.168 

0.053 

0.40 

8 

0.31  9.4 

71.2 

9.5 

3.30 

3.27 

323 

0.185 

0.068 

0.17 

11 

0.17  11.5 

6.48 

11.5 

1.40 

3.07 

333 

0.210 

0.090 

0.10 

14 

0.10  14.0 

1.41 

14.0 

0.660 

2.86 

343 

0.239 

0.135 

0.073 

24 

0.076  23.3 

0.272 

23.5 

0.240 

2.25 

353 

0.275 

0.192 

0.062 

34 

0.064  35.7 

0.139 

35.7 

0.130 

1.80 

363 

0.323 

0.255 

0.053 

45 

0.054  45.6 

0.097 

45.6 

0.095 

1.57 

373 

0.378 

0.310 

0.043 

47 

0.044  45.7 

0.079 

45.7 

0.078 

1.58 

Eg-High  Temp.  (kJ) 

18.6 

19.6 

43.6 

38.7 

Ejj-Low  Temp.  (kJ) 

63.0 

59.2 

257e 

68.5 

a  fit  by  Denault  and  Prud'homme  using  original  log-X^; 
b  recalculated  data  using  original  log-x^  and  non-integer  p’s; 
c  modified  log-x^,  equation  (4); 

<1  log-normal  distribution,  equation  (1 1); 

®  see  text  for  explanation. 


13 


O'Connor  and  Blum  L.og<x^  Distribution 

Figure  Caption 

Figure  1.  Xq  values  firom  the  log-normal  distribution's  fit  to  and  T2  data  (reference  (7))  as 
a  function  of  tempo'ature  and  PIPA  concentration.  The  curves  drawn  are  for  visual  clarification. 
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